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We discuss the strong-field dynamics of HeH + , the simplest stable heteronuclear molecule, focus- 
ing on identifying a laser regime for which there is a sufficient dissociation signal for experimental 
measurement. We numerically solve the time-dependent Schrodinger equation to obtain total dis- 
sociation probabilities, kinetic energy release spectra, and momentum distributions for wavelengths 
from 800 nm to 2400 nm. The suitability of this simple system as a prototype for understanding 
the strong- field nuclear dynamics of heteronuclear dissociation is discussed. 



I. INTRODUCTION 



There is considerable interest in the dynamics of 
molecules in intense laser fields [l], 0. Being the sim- 
plest molecules, and D^" have been at the center of 
many of these studies. Their simplicity allows for theoret- 
ical studies with minimal approximations [2- 9] , while the 
accuracy of these calculations opens the possibility for 
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quantitative comparison with experimental results 
Because of these features, and serve as bench- 
mark systems, providing insight into how more compli- 
cated molecules behave in intense laser fields. 

In these benchmark studies of and its isotopes, 
many interesting phenomena have been calculated and 
measured. Among these are above threshold dissocia- 
tion [ioj ]. below threshold dissociation 11 1, coherent con- 
trol of dissociation through variation of carrier envelope 
phase of a single color pulse [l2rlT4T | or the delay be- 
tween two pulses with different colors [l5l - [l7l |. and bond 
softening |18l - [2rj| |. Given the utility and success of these 
Hj" studies, a natural next step is to identify and study 
a similarly fundamental heteronuclear system Carrying 
out calculations on such a system will not only provide 
additional insight into how the aforementioned phenom- 
ena generalize, but it will also allow for the identifica- 
tion of new phenomena that are not present in homonu- 
clear molecules. We exclude the simplest heteronuclear 
molecule, HD + , because it does not produce an electronic 
asymmetry within the usual infinite nuclear mass Born- 
Oppenhcimer approximation. 

Thus, we propose HeH + as a benchmark system for 
studying the nuclear response of heteronuclear molecules 
to intense laser fields. The viability of this species as an 
experimental target has already been established. Specif- 
ically, ion beams of HcH + have been produced to demon- 
strate the existence of metastable HeH 2+ [2l|, [22j and, 
more recently, to serve as targets for an ultraviolet free- 
electron laser (23[. In this paper, we will focus on dis- 
sociation of HeH + at far- to mid-infrared wavelengths 
between 800 nm and 2400 nm. Intense, short pulses 
with such long wavelengths are becoming increasingly 
available experimentally and have been used in recent 
years to study high harmonic generation 24] and ioniza- 
tion [25L [2fj| . We will limit our exploration to five and ten 
cycle laser pulses with intensities between 10 12 W/cm 2 



and 10 14 W/cm 2 . Our results show that although the 
dissociation probability under these conditions is small 
for the more standard 800 nm pulses, the longer wave- 
lengths in our parameter space are capable of producing 
substantial dissociation, making this process experimen- 
tally accessible. 

We argue that the large difference in energy between 
the ground and first excited electronic channels for HeH + 
causes permanent dipole transitions to be much more 
probable than electronic excitation for the wavelengths 
considered. Because the system is dominated by sin- 
gle channel effects, it behaves analogously to an atom 
in a laser field. This atom-like behavior allows us to un- 
derstand physical observables calculated for HeH + using 
pictures previously developed to study intense field ion- 
ization. Moreover, it opens up interesting opportunities 
to observe phenomena usually associated with atoms in 
a molecular system. 



II. THEORETICAL METHODS 

A. Schrodinger Equation in the Single Channel 
Approximation 

All calculations in this paper use atomic units and are 
carried out by solving the time-dependent Schrodinger 
equation in the Born-Oppenheimer approximation. We 
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FIG. 1. Coordinates used in the Hamiltonian. 
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FIG. 2. (Color online.) (a) The adiabatic potential energy 
curves |27[ for the ground state of HeH + , X 1 E + , and the first 
excited singlet state, A 1 E + . (b) The X 1 E + potential used in 
our calculations [2g| plotted along with the vibrational lev- 
els for zero orbital angular momentum. The arrows indicate 
the energy after dissociation from the absorption of 800 nm, 
1200 nm, 1600 nm, 2000 nm, and 2400 nm (from left to right) 
photons from the vibrational ground state. 



treat the laser's electric field classically and use the length 
gauge within the dipole approximation. Using the dipole 
operator 
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the time-dependent Schrodinger equation is given by 
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where H e \ is the electronic part of the Hamiltonian in- 
cluding nuclear repulsion and £ is the electric field pro- 
duced by the incident laser pulse. In Eq. @, the coor- 
dinates are as shown in Fig. 1, and the reduced mass is 
given by 
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The Born-Oppenheimer potentials and electronic 
wavefunctions are found by solving 

ff c i$„ A (i?;rx ,r 2 ) = U vK {R) $„ A (Rpi ,r 2 ) . (4) 

In this paper, we use the potentials from Ref. [28j |. Be- 
cause <fr v \(R; r 1; r 2 ) form a complete set in the electronic 
coordinates, we can write 

*= ^F JaA (R,t)Dif A (4>,9,0)^ aA (R;r u r 2 ), (5) 
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where a denotes the electronic state, A is the magnitude 
of the projection of the electronic angular momentum on 
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FIG. 3. (Color online.) The body frame z-axis projection 
of the X X E+ permanent dipole, ($i|d|$i), and X 1 E+-4 1 E+ 
transition dipole, ($i|d|$2), from Ref. [23| along with our 
approximate permanent dipole (<3> |d. | <&) from Eq. j6]). 



the internuclear axis, J is the total orbital angular mo- 
mentum, M is the projection of the total orbital angular 
momentum along the lab frame z-axis, and D J MJ ^{(f>, 8, 0) 
is the Wigner £)-function normalized over the polar and 
azimuthal angles 9 and </> that describe the orientation of 
the molecular axis relative to the lab frame. 

In order to reduce the number of electronic channels 
required in Eq. (5), we examine the Born-Oppenheimer 
potentials and dipole couplings to understand the role 
of electronic excitation for our parameters. From Fig. [21 
we see that there is a difference of at least 0.40 a.u. in 
energy between the ground electronic channel X 1 Y. + and 
the first excited channel A 1 !]" 1 ". Therefore, a minimum of 
eight photons is required for electronic excitation with an 
800 nm pulse, the shortest wavelength that we consider. 
Moreover, we see from Fig. [3] that the magnitudes of the 
permanent dipole, (<I>i|d|<I>i), and the X 1 'E + -A 1 Y, + tran- 
sition dipole matrix element, (<I>i|d|$ 2 ), have nearly the 
same value for 1.0 a.u. < R < 2.0 a.u., and that the per- 
manent dipole is larger for R > 2.0 a.u. In the range 
where the matrix elements are comparable, however, a 
minimum of thirteen 800 nm photons is required for elec- 
tronic excitation. Taken together, these considerations 
suggest that electronic excitation will be negligible for 
all but very high intensities. Therefore, we will neglect 
electronic excitation in our calculations and keep only the 
state. Because our calculations only involve one 
electronic channel, we will drop the electronic channel 
labels for the remainder of this work. 

When our calculations were carried out, the extensive 
results for the dipole matrix elements in Ref. [27[ had not 
yet been published, so we used an approximate form of 
the permanent dipole: 
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This form is exact asymptotically as the system dissoci- 
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ates to He(ls 2 ) + H + and can be derived either classically 
or by direct evaluation in the separated atom limit. Fig- 
ure 3 compares this approximate dipole to the accurately 
calculated ones from Ref. [13] ■ The figure shows our ap- 
proximation to be in good agreement with the calculated 
dipole with discrepancies for small R. The fact that the 
actual dipole is smaller than the one used in our calcu- 
lations implies that the intensity required to produce a 
given dissociation probability will be larger than found 
in our calculations. However, lowest order perturbation 
theory predicts that the actual and reported intensities 
required to produce the lowest energy kinetic energy re- 
lease (KER) peak are of the same order of magnitude for 
the wavelengths considered. 

With the above approximations, Eq. (5) reduces to 
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where D^ Q — Yjq was used and we have assumed linearly 
polarized light in order to include only M = in the 
expansion. Substituting this sum into Eq. (2), projecting 
out $Yj7 , and using Eq. (6) gives the equations we solve 
for the radial part of the wavefunction: 



(8) 



where we have taken the laser polarization in the 
lab frame z direction. Equation (8) neglects all cou- 
pling between electronic states, including non-Born- 
Oppenheimer effects and Coriolis coupling |29| . 



B. Numerical Methods 




where to is the carrier frequency and r = tfwhm/v / 2Tti2. 
Our calculations are carried out for tfwhm correspond- 
ing to both five cycle and ten cycle pulses, using wave- 
lengths in the 800-2400 nm range and intensities between 
10 12 W/cm 2 and 10 14 W/cm 2 . Calculations begin at 
the time U when the intensity of the field is 10 7 W/cm 2 
and end at the time tf, after the peak intensity, when 
the intensity falls off to 10 s W/cm 2 . The time step is 
6 = 0.30 a.u., and we dynamically increase the number 
of partial waves included in our expansion of # 0, m, 
leading to a maximal expansion of between 10 and 27 
partial waves depending on the parameters of the laser 
pulse. For five cycle pulses, we use a non- uniform radial 
grid with i? m i n = 0.20 a.u., i? max = 120.0 a.u., and a grid 
spacing AR rj 0.00094 a.u. for 0.20 a.u. < R < 8.0 a.u. 
and AR w 0.0094 a.u. for R > 8.0 a.u. Three point 
averaging, where each of three adjacent grid points is 
weighted equally, is repeatedly applied 20,000 times to all 
of the grid points in order to smooth the abrupt change in 
grid spacing at i?=8.0 a.u. This smoothing is necessary 
for a non-uniform grid to avoid a loss of accuracy in the 
radial derivatives of our wavefunction near R = 8.0 a.u. 
For ten cycle pulses, we use the same grid spacing but 
extend R max to 200.0 a.u. Convergence with respect to 
all of these parameters was tested and found to give the 
KER spectrum accurate to three digits up to an energy 
of 0.20 a.u. 

Because part of our aim is to identify the parame- 
ters that give the largest dissociation signal, we choose 
the isotopes of helium and hydrogen that maximize the 
dipole moment and, thus, the dissociation at a given 
intensity. Equation (6) shows that in order to max- 
imize the dipole, we want the isotopes with the most 
asymmetric masses. Consequently, we consider 4 He and 
1 H. For their masses, we use vha = 7351.67 a.u. and 
ins = 1836.15 a.u. 



The remaining task is to solve Eq. (8). Time propa- 
gation of a given initial state is carried out according to 
Eq. (8) u sing the generalized finite differencing scheme 
from Ref. [301 ] to approximate the time-dependent Hamil- 
tonian as the block tridiagonal matrix H(i). Here, the 
diagonal blocks are the matrix representation of the field- 
free Hamiltonian, and each off-diagonal block is a diago- 
nal matrix with the J' to J' ± 1 dipole coupling. Using 
H (t) , time propagation is accomplished for small time 
steps 5 by 

F (R, t + 5) = e -im+s/2)8 F ^ t j _ (Q) 

where the components of F are the Fj in Eq. (7). The ac- 
tion of this short-time propagator is carried out by apply- 
ing the split-operator techniques and the Crank-Nicolson 
method outlined in Ref. [3lj . 

We use initial vibrational states with J — and express 
the electric field as 

£ = £ e- {t/T)2 cosut, (10) 



III. RESULTS AND DISCUSSION 

A. Dissociation Probability 

To quickly locate the parameters that produce sub- 
stantial dissociation, we first consider the total dissocia- 
tion probability 

^ = i-£l<x^(*/)>f. (ii) 

vJ 

for a series of initial vibrational states. In Eq. (11), Xvj 
are the field-free ro-vibrational nuclear eigenstates and 
Fj(tf) are the components given by Eq. (8). The in- 
tensity dependence of P for several initial vibrational 
states is presented in Fig. 4. From the figure, we see that 
P (x I n for most initial vibrational states, which is the 
result expected from lowest order perturbation theory for 
n photon absorption. Deviations from I n do appear for 
high intensities, indicating that lowest order perturbation 
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FIG. 4. (Color online.) Total dissociation probability P as a function of intensity for different wavelengths, pulse lengths, and 
initial vibrational states. The top figures are for 5 cycle pulses, while the bottom figures are for 10 cycle pulses. From left to 
right, the columns correspond to the initial vibrational states v=0, 3, 6, and 9. 
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FIG. 5. (Color online.) Bound-free dipole matrix elements, 
Eq. (12), as a function of KER E for v = 6 and v = 9. The 
horizontal arrows represent the range of energies expected for 
one photon absorption for wavelengths between 2400 nm and 
800 nm. 



theory is no longer adequate. Figure 2 shows that for an 
initial v = state, two photons are required for dissocia- 
tion at 800 nm and 1200 nm, three photons are required 
at 1600 nm and 2000 nm, and four photons are required 
at 2400 nm. These values of n are consistent with the 
slopes found in Fig. 4. This same analysis can be applied 
to explain the slopes seen for other initial states in the 
figure. 

From Fig. 4, we see that there do exist parameters 
that produce substantial dissociation: up to 8.5% (five 
cycles) and 16% (ten cycles) for 10 13 W/cm 2 , and up 
to 48% and 65% for 10 14 W/cm 2 , respectively. In all 
cases, these maximal P occurred for v = 6 at 2400 nm 
— not surprising considering only one photon is required 
for dissociation from this state. Dissociation from v = 9, 



however, also requires only a single photon, yet its disso- 
ciation probability is an order of magnitude smaller than 
for v = 6. The difference in dissociation probability for 
these two states is explained by their bound-free dipole 
matrix element, 
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where \E, J) is the the energy normalized field-free con- 
tinuum state with energy E and angular momentum J. 
Figure [5] shows De9 and Deg together with the ranges of 
KER expected in each case for the one-photon peak for 
wavelengths between 2400 nm and 800 nm. Keeping in 
mind that the vertical scale in the figure is logarithmic, 
it is clear that De9 is considerably smaller than Des at 
any given wavelength, yielding, in turn, a much smaller 
dissociation probability. 

The wavelength dependence of P evident in Fig. @] can 
also be explained by Dev in Fig- we see that for 
v=6 and 9 the dipole matrix elements decrease mono- 
tonically with increasing KER. In fact, this behavior is a 
general feature of De v for this system. Therefore, longer 
wavelengths are expected to give the largest dissociation 
probabilities since they produce dissociating fragments 
with the lowest KER. 

There is also substantial dissociation for other initial 
states and laser parameters. Both five and ten cycle 
pulses with / = 10 14 W/cm 2 produce dissociation prob- 
abilities greater than 1.0% for both A = 2000 nm and 
A = 2400 nm in all studied initial vibrational states ex- 
cept v = 0. Because the largest dissociation occurs for 
intensities between 10 13 W/cm 2 and 10 14 W/cm 2 and 
wavelengths between 2000 nm and 2400 nm, we constrain 
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our studies of the KER and momentum distribution to 
this region of parameter space. 



B. Kinetic Energy Release Spectrum 

Due to the fact that HcH + is a benchmark molecular 
system, it is useful to try to understand its behavior in 
a laser field using the pictures already developed for the 
benchmark molecule H^". In particular, the dressed po- 
tentials of the Floquet representation have proven espe- 
cially useful in understanding the kinetic energy release 
spectrum of H^" [l|, H, HH . These dressed potentials are 
obtained by shifting the Born-Oppenheimer potentials 
by integer multiples of u>, keeping only those states con- 
nected by one or more application of the dipole selection 
rules. At crossings of the dressed potentials, transitions 
are much more likely because these crossings correspond 
to a resonance-like condition. 

Because there is only one electronic channel relevant 
for our laser conditions, shifting it by multiples of uj pro- 
duces parallel potentials with no crossings. When the 
appropriate centrifugal barriers are included, there are 
crossings, but only for such high J and photon number 
that they are not likely to produce transitions. Thus, 
non-resonant transitions must dominate for the laser pa- 
rameters we consider, reducing the utility of the dressed 
potential picture in this study. We note that this fail- 
ure of the dressed photon picture will be generally true 
for processes dominated by permanent dipole transitions 
and is not a feature specific to HcH + . If we were con- 
sidering shorter wavelengths, such that transitions to the 
A 1 Yi + (or higher) state were relevant, then the dressed 
potential picture would again prove useful. 

It seems, then, that to understand the dissociation of 
HeH + , we must think about the process differently than 
for H^. A natural picture to adopt is that for atomic 
ionization because including only a single channel makes 
the HeH + nuclear time-dependent Schrodinger equation, 
Eq. (8), look just like that of a hydrogenic atom. Pur- 
suing this analogy, Eq. (8) shows that the "electron" in 
HeH + has a charge of 171a/ (niA + tub), a mass of hab, 
and interacts via a short-range central potential. This 
analogy allows us to interpret our results using the pic- 
tures that are used to understand intense field ionization. 

Atomic ionization is primarily understood using one 
of two pictures: tunneling or multiphoton ioniza- 
tion. Which picture is most appropriate is customar- 
ily determined by the value of the Keldysh parameter 
7 = \J Eb/2U P , where Eb is the binding energy and 
Up~q • E 2 1 'Afiuj 2 is the ponderomotive energy [32|,[33j|. In 
the case of ionization, the charge q and mass fi in U p 
refer to an electron and are unity in atomic units; they 
will differ substantially from unity, though, for our analog 
system. For HeH + (v=3, J=0) in a 2400 nm laser pulse 
with peak intensity 7=8.9xl0 13 W/cm 2 , 7=4.6, plac- 
ing the system in the multiphoton regime. For higher- 
lying ro-vibrational states, the Keldysh parameter will 
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FIG. 6. (Color online.) The KER spectrum for an initial 
vibrational state of v — 3 in a five cycle pulse The vertical 
dashed line at E/ui — 0.26 corresponds to the expected loca- 
tion of the first peak, and the dashed diagonal line shows the 
dP/dE oc I n behavior expected from lowest-order perturba- 
tion theory. 



become smaller at these laser parameters, dropping to 
7=0.67 for HeH + (t>=9, J=0). So, with relatively modest 
— and likely achievable — changes in the laser parame- 
ters, HeH + can studied over the whole range of phys- 
ical regimes from multiphoton to tunneling. For the 
remainder of this paper, however, we will restrict our 
consideration to HeH + (i>=3, J=0) in a 2400 nm pulse 
with 7=8.9xl0 13 W/cm 2 which lies in the multiphoton 
regime. Consequently, we expect to see many of the same 
phenomena that are seen in the multiphoton ionization 
of atoms, such as above threshold ionization (ATI), in 
the strong-field dissociation of HeH+ . 

One commonly studied observable in atomic systems 
that reveals these phenomena is the photoelectron spec- 
trum. In HeH + , the analog is the KER spectrum. That 
is, the probability of the molecule dissociating with a rel- 
ative nuclear kinetic energy E, 

dP 

— = J2\(E,J\Fj(t f ))\ 2 . (13) 
J 

A typical result for the KER spectrum in the parameter 
space that we are considering is shown in Fig. 6. This 
spectrum clearly shows the characteristic photon-spaced 
peaks expected for above threshold ionization (ATI) - 
except, of course, that this is the nuclear KER spectrum, 
making it above threshold dissociation (ATD). In the 
atomic case, the energies of the ATI peaks are given by 
E n = nu — Eb — Up, where n is the net number of photons 
absorbed. Because the ponderomotive energy is negligi- 
ble for the parameters we are considering (U p /w = 0.040 
for the parameters in Fig. 6 and Fig. 7), this result sim- 
plifies to E n /bj = Eb/cj + n. For v = 3 and A = 2400 nm, 
we expect the first peak to occur for two photon absorp- 
tion at E n= 2/u) = 0.26, which is in good agreement with 
the result in Fig. 6. 

It is also important to note the large number of ATD 



peaks visible in Fig. 6. Our result shows more ATD or- 
ders than is typical and is a consequence of the large 
energy difference between the A 1 !]"*" and poten- 
tials, which allows for the absorption of many photons 
before electronic excitation plays a role and blurs the 
peaks. Even so, the highest ATD peak shown in Fig. 6 
lies 0.19 a.u. — more than eleven photons — below the 
minimum of the A 1 T, + potential. 

One notable difference between Fig. 6 and a typical 
ATI spectrum is the absence of the plateau usually seen 
for 3.17U p < E < 10U P (32)]. Because the ponderomo- 
tive energy is small relative to the frequency of the laser, 
our ATD plateau would occur between E/ui = 0.13 and 
E/ui = 0.40. The fact that lQU p /u - 3.17U p /uj is much 
smaller than one explains our inability to observe an 
ATD plateau. Physically, this small interval is due to 
the fact that the massive nuclei are unable to gain a sub- 
stantial amount of energy in the field. For a pulse with 
/ = 8.9 x 10 13 W/cm 2 , it would require a wavelength of 
almost 3.7 /im to produce a plateau region with a large 
enough energy range to fit two KER peaks. 

Finally, the spectrum in Fig. 6 shows that high-order, 
non-linear processes are occurring: the figure gives evi- 
dence for at least eleven photon absorption. Neverthe- 
less, Fig. 4(b) suggests that the deviation from lowest- 
order perturbation theory is small. This apparent in- 
consistency can be reconciled by recognizing that the 
area under the n = 2 peak dominates the energy integral 
needed to calculate dissociation probability from dP/dE. 
However, Fig. 6 also shows that the non-perturbative 
effects extend beyond the presence of high-order ATD 
peaks. The diagonal line indicates I n behavior that 
matches the highest energy peaks as expected. The lower 
peaks, however, do not meet this line, suggesting the con- 
tribution of high-order pathways in these peaks. Due to 
their dominating the energy integral of dP/dE, it is this 
deviation from I n that produces the non-perturbative be- 
havior seen in Fig. 4(b) rather than the presence of the 
higher-order peaks. 

C. Momentum Distribution 

We can also use the atomic analogy to understand the 
distribution of relative nuclear momenta k 

^ r ^m f ))\\ (14) 

where Ok is the scattering angle with respect to the laser 
polarization and the state |k) is an energy- normalized 
scattering state with asymptotic out-going momentum k. 
Note that Eq. (fl3|) is obtained from Eq. ([M)) by integrat- 
ing over the angle Ok- The momentum distribution for 
v = 3 and A = 2400 nm, the same case shown in Fig. 6, is 
shown in Fig. 7. We, of course, still see the photon-spaced 
peaks with decreasing probability, as expected from the 
atomic ionization picture and Fig. 6. The angular distri- 
bution of each ATD peak is consistent with the identifica- 
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FIG. 7. (Color online.) "Momentum" distribution from 
Eq. (|14|l plotted as a function of KER E and angle 9k for 
an initial ro-vibrational state with v=3 and J=0 in a five 
cycle pulse with an intensity of 8.9xl0 13 W/cm 2 and a wave- 
length of 2400 nm. The upper-axis label n indicates the net 
number of photons absorbed. 



tion of the number of photons involved. For instance, the 
peaks at E/u) = 0.3,2.3,... are nonzero at 9k/ir = 0.50 
as expected for even parity final states produced by the 
absorption of n = 2,4, ... photons from the J = initial 
state. Similarly, the peaks at E/u = 1.3,3.3,... have a 
node at 9k /V = 0.50 corresponding to an odd number of 
photons absorbed — specifically n = 3, 5, ... — from the 
J = initial state. 

In both Fig. [5] and Fig. there is additional struc- 
ture on the two-photon peak in the form of a peak at 
E/lo = 0.50 (indicated by arrows in the figures). This 
peak can be explained by noting that the v = 3 to 
v = 6 transition is nearly on resonance for 2400 nm: 
W63=1.2o;. Even though the effective intensity at CJ63 
is reduced by seven orders of magnitude, the resonant 
enhancement is still sufficient to produce a peak. More- 
over, this resonance enhanced multiphoton dissociation 
(REMPD) mechanism correctly predicts the peak to be 
at E/ui = 0.50 in agreement with our calculations. A 
REMPD peak is noticeable on the n = 3 peaks in the 
figures as well. For higher order peaks, however, the 
REMPD peak grows increasingly broad in energy due 
to the bandwidth of the laser pulse, washing it out as 
a distinct peak. As its name suggests, REMPD is the 
analog of the well-known resonance enhanced multipho- 
ton ionization in ATI spectra [33j , and we expect it to be 
ubiquitous in long wavelength dissociation of HeH + since 
these wavelengths are more likely to drive a bound-bound 
resonance. 



IV. SUMMARY AND OUTLOOK 

We have proposed HeH + as a benchmark system for 
understanding the nuclear dynamics of heteronuclear 



7 



molecules in intense laser fields. However, given that 
dissociation of HcH + at the common intense laser wave- 
length of 800 nm is consistently below useful experi- 
mental detection levels at intensities that will not also 
ionize the molecule, it was necessary to explore other 
laser parameters. We find that increasing the wave- 
length is an effective way to increase the dissociation 
probability at modest intensities. In particular, for wave- 
lengths between 2000 nm and 2400 nm and intensities be- 
tween 10 13 W/cm 2 and 10 14 W/cm 2 , we see probabilities 
greater than one percent for v = 3, 6, and 9, making dis- 
sociation of HeH + experimentally accessible. In general, 
we expect even larger HeH + dissociation probabilities for 
wavelengths longer than the 2400 nm considered here. 

It turns out that HeH + is a rather unique molecule 
in that the minimum electronic excitation energy is rel- 
atively large. While this property decreases its utility as 
a prototype for the behavior of heteronuclear molecules 
generally, it does allow for the surprising possibility of 
understanding the nuclear dynamics of HcH + in terms 
of atomic ionization pictures. The same reduction of the 
system to a single channel also has the intriguing conse- 
quence that all of the nuclear dynamics for a substantial 
range of laser parameters is induced by the molecule's 
permanent dipole. This lies in contrast to the electronic 
transition dominated strong-field response of most every 
other molecule that has been studied to date. 

So, while HeH + may not be the general prototype 
we aimed for, it is still a very interesting molecule to 
study. Even though much of the physics for this system 
can be understood using atomic pictures, HeH + still ex- 



hibits molecular behavior that differentiates it from the 
atomic systems that have been studied in the past. As a 
molecule, for example, the initial state of the system will 
actually be an incoherent distribution of ro-vibrational 
states determined by the temperature and mechanism for 
creation of the molecule — unlike atoms. The strong- field 
dynamics of HeH + is thus not exactly molecule-like, and 
it is not exactly atom-like. But, with only two electrons, 
it is a benchmark system for which precision calculations 
and experiments ought to be possible to elucidate the 
dynamics of a system that displays physical behavior in 
limbo between that of molecule and that of an atom. 
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